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Intermediate regimes in granular Brownian motion: Superdiffusion and subdiffusion 



(N 

o 

ov 



c3 
c3 



o 
o 



>: 

VO, 
OV 

o: 



X 



Anna Bodrova, 1 Awadhesh Kumar Dubey, 2 Sanjay Puri, 2 and Nikolai Brilliantov 3 

'Faculty of Physics, M.V.Lomonosov Moscow State University, Moscow, 119991, Russia 
2 School of Physical Sciences, Jawaharlal Nehru University, New Delhi, 110067, India. 
'Department of Mathematics, University of Leicester, Leicester LEI 7RH, United Kingdom 

(Dated: September 20, 2012) 

Brownian motion in a granular gas in a homogeneous cooling state is studied theoretically and 
by means of molecular dynamics. We use the simplest first-principle model for the impact-velocity 
dependent restitution coefficient, as it follows for the model of viscoelastic spheres. We reveal that for 
a wide range of initial conditions the ratio of granular temperatures of Brownian and bath particles 
demonstrates complicated non-monotonous behavior, which results in transition between different 
regimes of Brownian dynamics: It starts from the ballistic motion, switches later to superballistic one 
and turns at still later times into subdiffusion; eventually normal diffusion is achieved. Our theory 
agrees very well with the MD results, although extreme computational costs prevented to detect 
the final diffusion regime. Qualitatively, the reported intermediate diffusion regimes are generic for 
granular gases with any realistic dependence of the restitution coefficient on the impact velocity. 
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Introduction. Brownian motion is a fundamental pro- 
cess in nature which takes place on very different time 
and length scales In its classical formulation it im- 
plies a random motion of a big (Brownian) particle driven 
by thermal motion of much smaller bath particles. The 
mean-square displacement of Brownian particles grows 
with time as (R%) ~ t 13 , with /3 = 1 for normal diffu- 
sion. Anomalous diffusion with /3 > 1 (superdiffusion) 
and j3 < 1 (subdiffusion) has been found in a variety of 
systems 0, ranging from supercooled and glass forming 
liquids, e.g. [3|, to surfaces [J], biological [5( and gran- 
ular [6[ systems; usually a few diffusion regimes with a 
crossover between them are observed. Mechanisms of 
anomalous diffusion are seemingly very different - macro- 
molecular crowding in solutions and cells caging in 
glassy liquids [3J], trapping in a random potential [J], etc. 
Nevertheless, one can possibly infer an important com- 
mon feature - anomalous diffusion always arises due to 
interplay of different basic dynamics in a system: For in- 
stance, a (fast) hopping dynamics in glass forming liquids 
interferes with a (slow) dynamics of structural rearrange- 
ments, resulting in subdiffusion. Moreover, anomalous 
diffusion in most of the systems manifests as an interme- 
diate regime, which transforms asymptotically into nor- 
mal diffusion, e.g. [HI]. 

In the present Letter we report intermediate anoma- 
lous diffusion in granular systems, where particles inter- 
act with dissipative forces. Although the microscopic 
mechanism of this phenomenon is novel and related to 
inelastic collisions, it again results, as it is shown below, 
from the interplay of different dynamics represented by 
evolution of two granular temperatures. 

Granular Brownian motion demonstrates a large va- 
riety of new surprising phenomena, among which is the 
violation of energy equipartition in a mixture of mas- 
sive and light particles [7H9|: Both granular tempera- 



tures of Brownian (Tb) and bath (T) particles decrease 
in a force-free granular gas, while their ratio rapidly re- 
laxes to a steady-state value (Tb/T) s . s . > 1. Moreover, 
(Tb/T%. s . increases with increasing inelasticity and the 
mass ratio m^/m of Brownian (tub) and bath (m) parti- 
cles @, IH • It has been also shown that Brownian motion 
in these systems may change qualitatively - from diffu- 
sive to a ballistic one, depending on the steady state value 
of (p = (Tgm) / (Tms) [10) . Surprisingly, one can treat 
this phenomena as a "phase transition" with the "order 
parameter" ip [Io| . which however is not observed as a 
real process, since the steady state (p is constant. 

To date, the most theoretical studies of Brownian mo- 
tion are done with the assumption of a constant resti- 
tution coefficient e = g//g Q], where g and gf are, re- 
spectively, the normal components of the relative inter- 
particle velocities before and after a collision. The sim- 
plifying hypothesis of constant e, however, contradicts 
experimental results, which indicate that e does depend 
on g and that e(g) is a decreasing function of the impact 
velocity ; it is not also consistent with a rigorous the- 
oretical analysis 0- The impact- velocity dependence of 
the restitution coefficient may, however, crucially change 
the behavior of a granular system, e.g. [121]. In particular, 
it has been shown that in a gas of viscoelastic particles, 
where the impact- velocity dependence of e stems from the 
simplest first-principle model 13 [jClusters and vortices 
manifest as transient phenomena [l4| , in a sharp contrast 
with the case of a constant e 15|. Furthermore, in a mix- 
ture of massive and light particles, the temperature ratio 
Tb/T, as well as (p, evolve with time in a HCS [la ], again 
in a contrast with the case of e = const. Therefore, one 
can expect that in a granular gas of viscoelastic particles, 
where ip varies with time, the "phase transition" between 
different regimes of Brownian motion may occur during 
the gas evolution as a real physical process. 
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We analyze granular Brownian motion in a gas of vis- 
coelastic particles in a homogeneous cooling state (HCS) 
- the basic state of a force free granular gas by means of 
molecular dynamics (MD) and theoretically. We have re- 
vealed a sequence of new diffusive regimes for a Brownian 
particle - superballistic motion (superdiffusion), which 
follows after ballistic one and transforms then to subdif- 
fusion. Eventually, normal diffusion is observed. 

MD simulations. We perform event-driven simula- 
tions [I?} for a force-free system, using N = 64000 bath 
particles of mass m and diameter a and one Brownian 
particle of mass tub 3> m and diameter o~b (for sim- 
plicity we take erg = a) in a cube of length L = 130 a 
with periodic boundary conditions. The number density 
of the bath particles was n = TVttct 3 /(6i 3 ) ~ 0.015. The 
reported results correspond to averages over 20 or 50 in- 
dependent runs, depending on the initial conditions. 

In our simulations we use the restitution coefficient^ as 
it follows for the model of viscoelastic particles IH 1 
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Here, Cfe's are numerical coefficients, which have been 
computed up to k — 20 [HI and 5 quantifies the dissi- 
pative interactions. Eq. (JlJ describes both types of col- 
lisions - between Brownian particle and a bath parti- 
cle and between the bath particles themselves. In the 
former case w = \(cb% ■ e)|, where cbi = VBi/vr, with 
VBi = vb ~ Vi being the relative velocity of Brownian 
(vb) and i-th bath (v*i) particle, vt — y^2T/m is the 
thermal velocity of the bath particles and e is the unit 
vector, joining particles' centers at the collision instant. 
In the latter case w = \(c\j ■ e)\, where Cij = Vij/vT, with 
Vij = Hi — Vj being the relative velocity of the colliding 
i-th and j-th bath particles. For both types of collisions 
we use for simplicity the same 6 = 0.1. The quantity 
u(t) = T(t)/T(0) is the dimensionless temperature of the 
bath at time t, and granular temperatures are defined as 
usual, 3T/2 = (mv 2 /2) and 3T B /2 = (m B u|/2) . The 
results of the MD simulation are presented in Fig. [T] 

Theory. Since concentration of Brownian particles is 
much smaller than that of the bath particles, we as- 
sume, that Brownian particles do not affect evolution of 
the bath Q: In a force- free gas T gradually decreases, 
due to dissipative collisions, with the cooling coefficient 
f (t) = - (dT/dt) /T (t). As temperature T tends to zero, 
e tends to unity, and all collisions tend to be elastic 0- 

Evolution of Brownian particles may be described, us- 
ing the Boltzmann equation, /s(us,t) = /(/bJ) for 
the velocity distribution function, /s (??b , i) , e.g. [7|; the 
collision integral 7(/g,/) accounts for the alteration of 
fB(vs,t) in pair- wise collisions with the collision rule, 



vb - (1 + e)/j, (vsi ■ e) e. 



(2) 



collision and a = — —, — . Solving the Boltzmann equa- 

r" m B +m t> l 

tion with the standard technique, e.g. [7|, one obtains 
the distribution function fg. Brownian dynamics may 
be described with very compact notations if the pseudo- 
Liouville operator L is used 0: 

v B = Lv B = L v B + TbiVb , (3) 

i 

where Lq — vb^t b describes the free streaming and 

T B i = OQ JdeO (- (v B i ■ e)) \v Bl • ej 5 (rBi~ea ) (b Bl - 1) 

the binary collisions of particles. The radius- vector fsi 
in the above equation joins centers of Brownian and bath 
particle, Co = (ctb + o) /2 and the operator bsi acts on 
the velocities, as bsiVs = v' B with the collision rule ([2]). 

Taking into account, that i?s(t) = v B (ti)db\, 
the mean-square displacement of Brownian parti- 
cles reads in terms of the velocity autocorrelation 
function: = j* d*i J dt 2 (v B (h) v B (h)). 

We introduce the reduced time tb, defined as 
dr B = dty/ T B {t)/T B {fi)lT c {Q), where r" 1 ^) = 
A^/TTa 2 g2(o')n^jT{t)/m is the mean collision time of bath 
particles, with 52(c) being the contact value of the pair 
correlation function 7]. In this time scale one deals 
with the reduced velocities, Cb = v B / y/ 2Tg /m B , and 
granular temperature of Brownian particles remains con- 
stant. Hence the decay of the (reduced) velocity corre- 
lation function occurs similarly, as in equilibrium gases: 
A particle loses memory of its initial velocity in random 
collisions; in dilute gases the velocity correlation function 
is exponential [l9| . therefore we approximate 

(c B Mil)) cb Mfc))) = (4) exp f TB[h) ~ TB[h) 



T V B (h) 



(4) 

where (c%) =3/2 and f vB is the (reduced) velocity relax- 
ation time. Exploiting the above form of the correlation 
function, we obtain the mean-square displacement: 



=6 / dhD B {ti) 



1 — exp 



TB(t) - T B {tl) 



TvB (*l) 

Here Db = T B (t)T vB (t) /tub is the time-dependent diffu- 
sion coefficient of a Brownian particle, expressed in terms 
of T vB (t) = f vB (t)T c (0)/y/T B (t)/T B (0) - the velocity re- 
laxation time in laboratory time units. The reduced re- 
laxation time t v b (t) in the exponential function in Eq. (|4]) 
may be easily found from its time derivative at t\ — > 
yielding @, 



1 ( > (vbT Bl Vb) I, 
r vB =-(N-l) ( . 2) 

with the cooling coefficient of Brownian particles 



Here v' B is the velocity of a Brownian particle after a 
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Performing the averaging in Eqs. (0) and ©, we ignore, 
for simplicity, deviations of the velocity distribution func- 
tion from the Maxwellian 12[ ; it may be however shown 
that their impact on the calculated quantities is negligi- 
ble. Hence, we obtain the velocity relaxation time 



r~s(t) [8u~fi 2 (l + <p) 3/2 cr 2 g 2 (o- ) 
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and the cooling coefficient of Brownian particles: 



fc(t) = 2r- 1 (0)\/-^(l + ^) 1 



1- A* 
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<r 2 g 2 (a) 
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where A = 4C; + 53j+/»=i Q/Cfe are P ure numbers and 

Bi {t) = #/* (2u(t)f™ (i+mr /2 ° ( 2 -^pj r (| 



with r (x) being the Gamma-function. 

The temperature of the Brownian particles Tb can 
be found from the equation dTB(t)/dt — — £s(i)lB(i), 
while the temperature T(t) and the mean-squared dis- 
placement (i? 2 (i)) of the bath particles - from the above 
results using the substitute ms — >• m and <tb — > c. 

Results and discussion. The theoretical predictions 
are compared with the simulation data in Fig. [TJ First 
we consider the case of the initial energy equipartition, 
T B (0) = T(0) = 400/3. The time dependence of the ratio 
Tb /T demonstrates here a complicated non- monotonous 
behavior: It initially increases, reaches a maximum, and 
then decreases, tending eventually to unity, Fig.QJ,. This 
is in a sharp contrast with the case of a constant e, where 
Tb/T rapidly relaxes to a steady-state value. One can 
explain this effect as follows. From the collision rules 
[see, e.g. Eq. |(5J|], the ratio of energy losses of Brown- 
ian (AEb) and bath (AE) particles at a collision, scales 
as AEb/AE ~ vb/v, ~ y/ (T B m)/(Tm B ). Initially 
AEb -C AE, since Tb = T and raj 3> m, therefore 
Brownian particles cool down significantly slower than 
the bath particles and the temperature ratio Tb /T in- 
creases. One can call this effect as a "retarded cool- 
ing", which may be also understood as a tendency of 
a system to reach the ratio of Tb/T > 1, correspond- 
ing to the steady-state value (Tb/T) s . s . > 1 of a gas 
with the respective constant e. The larger the mass ratio 
the more pronounced the effect (Fig. [T^) . In the course 
of time the temperature ratio Tb/T gets so large that 
AEb and AE become comparable; the granular temper- 
atures start then to equilibrate. At this stage Tb/T de- 
creases - Brownian particles cool faster than the bath 
particles and one can call this regime as "accelerated 
cooling". Again, this effect may be understood, notic- 
ing that the ratio of (Tb/T) s s . > 1 is smaller for larger 
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FIG. 1: (a) Time dependence of the temperature ratio of 
Brownian (Tb) and bath (T) particles, (b) The mean-square 
displacement of Brownian particles (R%) as the function of 
the re-scaled time r u and the laboratory time t (inset). Lines 
- theory, symbols - MD data. Numbers indicate different 
regimes: 1 - ballistic motion, (R%) ~ T u\ 2 - superballistic 
regime (superdiffusion), (R%) ~ with /3 > 1; 3 - subdiffu- 
sion, (R%) ~ t„ with /3 < 1; 4 - normal diffusion (R%) ~ T u- 



e and that the effective restitution coefficient grows as 
the gas cools down, see Eq. ([1]). For very large times e 
tends to unity, the system becomes elastic and the energy 
equipartition is achieved. This complicated temperature 
dependence results in transition between different diffu- 
sion regimes, which may be most clearly seen using the 
re-scaled time r u , defined as dr u — dty/u(t) /t c (Q). In 
this time scale the average velocity of the bath particles 
is constant and they move, as in an equilibrium molecular 
gas, that is, ballistically, (R 2 ) ~ t 2 , for t u ~ 1 and dif- 
fusively, (R 2 ) ~ T„, for t u ^> 1 Brownian motion 
demonstrates, however, four different regimes - ballistic 
(regime 1 in Fig.[T|D), superballistic motion or superdiffu- 
sion (regime 2), subdiffusion (regime 3) and normal dif- 
fusion (regime 4). Physically, the superballistic motion, 
which starts after a short ballistic regime corresponds to 
the "retarded cooling", while the subsequent subdiffu- 
sion - to the "accelerated cooling". In the former case 
Brownian particles get hotter and hotter with respect to 
the surrounding gas and (R'b) grows with time faster 
than for normal diffusion; in the latter case, Brownian 
particles cool more rapidly than the gas, so that (R%) 
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grows with time slower. The duration of the subdiffusion 
regime increases with the mass ratio ms/m and may 
be very long. Asymptotically, at r„ — > oo, the system 
returns to the equipartition, demonstrating the normal 
diffusion, (R%) ~ t u . 

The respective transition between diffusion regimes is 
also seen for the laboratory time t (Fig.QJ), inset), where 
the "normal diffusion" is described by (R%) ~ t 1 ^ 6 , as 
it follows for the gas of viscoelastic particles @, 0] • If 
we consider the time scale tb, where Tb keeps constant, 
the superballistic motion does not appear; all three other 
diffusion regimes, nevertheless, present. 

Since superballistic motion is caused by the "retarded 
cooling", it is not observed if Tb/T does not increase. 
This happens for initial conditions, when the ratio Tb/T 
is already large, or more precisely, when Tb/T exceeds 
some threshold, estimated as the ratio (Tb/T) ss for a 
gas with a constant e, equal to the effective restitution 
coefficient e e g for a gas of viscoelastic particles at t = 0. 
In our case Eq. (TT]) with u(0) = I and (w) — 4/v2~7T yields 
£ eff (0) = 0.876 and T s (0)/T(0) = {T B /T) S . S . = 213.5. 
Hence, this initial temperature ratio delimits two incip- 
ient diffusion regimes: For T B (0)/T(0) < 213.5 the bal- 
listic and then superballistic motion takes place, while 
otherwise only ballistic regime is observed, Fig. QJ). In- 
terestingly, the ballistic motion in the latter case lasts for 
significantly longer time, than for normal diffusion; this 
follows from the fact that trajectories of heavy Brownian 
particles are almost not affected by much colder bath of 
light particles. The subsequent subdiffusion and normal 
diffusion persist independently on initial conditions, see 
Fig. Hb. This is the universal feature of Brownian mo- 
tion in a gas of viscoelastic particles. As it follows from 
Fig. [U a very good agreement between our theory and 
simulation data is observed for evolution of temperature 
and mean-square displacement for intermediate ballistic, 
superballistic and the onset of subdiffusive regimes; due 
to extreme computational costs we fail to detect numer- 
ically the final regime of normal diffusion. 

Conclusion. We have studied Brownian motion in a 
force- free granular gas, composed of particles interact- 
ing with an impact-velocity dependent restitution coef- 
ficient s(g), as it follows from the model of viscoelas- 
tic spheres. We have revealed that for a wide range of 
initial conditions the ratio of granular temperatures of 
massive and light particles Tb/T demonstrates compli- 
cated non-monotonous dependence on time as the sys- 
tem evolves in the homogeneous cooling state. Some- 
what similar to other systems with anomalous diffusion, 
this interplay of different dynamics of Brownian and bath 
particles, quantified by their temperatures, gives rise to 
a sequence of intermediate diffusive regimes in granular 
Brownian motion: At early times the ballistic motion 
is observed, which alters then to superballistic one (for 
non-monotonous evolution of Tb /T) . At still later times 
the opposite regime of subdiffusion - motion, slower than 



normal diffusion, always takes place; it lasts for a rela- 
tively long time and tends asymptotically to normal diffu- 
sion. Qualitatively, these features of granular Brownian 
motion are not specific to a particular model of e{g) but 
generic for granular systems with any realistic restitution 
coefficient that increases with decreasing impact velocity. 
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